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Abstract 



We study the decay B c — > D*j. There are two mechanisms contributing to the 
process. One proceeds through the short distance b^s^f transition and the other 
occurs through weak annihilation accompanied by a photon emission. The electro- 
magnetic penguin contribution is estimated by perturbative QCD and found to be 
4.68 x 10~ ls GeV . In particular, we find the contribution of the weak annihilation is 
6.25 x 10 _18 GeV^ which is in the same order as that of the electromagnetic penguin. 
The total decay rate T(B C — ► D*j) is predicted to be 1.45 x W~ 17 GeV and the 
branching ratio Br(B c — ► D*^f) is predicted to be 2.98 x 10~ 5 for tb c = 1.35ps. The 
the decays B c — > D*7 can be well studied at LBC in the near future. 
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1 INTRODUCTION 



Recently the physics of B c meson have got intensive attention [1,2]. It is believed that the 
B c meson is the next and the final member of B mesons. It provides unique opportunities 
to examine various heavy quark fragmentation models, heavy quark spin-flavor symmetry, 
different quarkonium bound state models, and properties of inclusive and exclusive decay 
channels. Being made of two heavy quarks of different flavors, B c weak decays also offer 
a rich source to measure CKM matrix elements of the standard model [2]. 

The interest of studying rare B decays lies on the fact that those decays, induced 
by the flavor-changing neutral currents b — > sj, are usually controlled by the one-loop 
electromagnetic penguin diagrams. They play important role in testing loop effects in the 
standard model and in searching for physics beyond the standard model( so called new 
physics). Most recently the weak radiative decays of B meson and bottom baryons have 
been systematically studied in [3,4]. Hence a precise study of the weak radiative decays 
of B c needs to be undertaken. In the present paper, we will address B c radiative decay 
B c -> D* sl . 

The subprocess b — > sj, taken as a free decay, is usually treated as the only flavor- 
changing contribution to B — > if*7[5], so do B c — > D*j. However, boundstate effects 
could seriously modify the results of the assumption. Boundstate effects include modifi- 
cations from weak annihilation which involves no flavor-changing neutral currents at all. 
The effects of the weak annihilation mechanism are expected to be large in the B c decays. 
We shall address this point in detail below. 

Unfortunately, the well known chiral-symmetry[6] and the heavy quark symmetry[7] 
cannot be applied straightforwardly to this process. Recently, a perturbative QCD(PQCD) 
analysis of the B meson exclusive decays seems give a good prediction [8]. As it is argued 
in ref. [9] that B c two body nonleptonic decay can be conveniently studied within the 
framework of PQCD suggested by Brodsky-Lepage[10] and then developed in [8]. 

The reason to use PQCD to analyze of B c radiative decay is as bellow. In the process 
b — > S7, s quark obtains large momentum by recoiling. In order to form a bound state 
with the spectator c quark, the most part of the momentum of s quark must be transferred 
to c by a hard scattering process, since in the final bound state (i.e D*~) the heavy charm 
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should share the most part of the momentum of D*~ . The hard scattering is suitable for 
PQCD calculation [8, 10]. The Feynmen diagrams we will calculate are given in Fig.l. 

The paper is organized as follows. In section 2, we display our calculations. We present 
our numerical results in Section 3. Section 4 contains the discussions and conclusions. 



2 CALCULATION 

Exclusive processes at large momentum transfer are exploited by Brodsky-Lepage[10] 
within PQCD starting with a Fock component expansion of the involved hadrons, where 
a twist expansion suggests that the contribution from the lowest order Fock component 
dominates the physical observable under consideration. An exclusive process then involves 
a perturbatively calculable hard scattering amplitude convoluted with nonperturbative 
soft physics wavefunctions of the initial and the final hadrons. Although these wavefunc- 
tions are incalculable from first principle, but they are universal for each hadron, i. e they 
are factorized out from the hard scattering amplitude and hence are independent of the 
process involved. 

The factorization scheme advocated by Brodsky-Lepage[10] is employed, where the 
momenta of quarks are taken as some fractions x of the total momentum of the meson 
weighted by a soft physical distribution functions 4>h{x). The peaking approximation is 
used for 0^(x)[ll]. The distribution amplitude of B c is 

0Bc(z) = ^y=fB c S(xi ^ ), (1) 

2V3 mi + m 2 

where rrii corresponds to m&, m c . The distribution function of D* is 

^(.H-i^-ig-), (2) 



where e^* is defined as 



e D * s = 2 • 3 



The spinor part of B c and the final meson D* are 



** +m <S 5 , (4) 



V2 -J2 
where e is the polarization vector of D*. 



2.1 Electromagnetic penguin contributions 

Within the standard model, the process is governed by the electromagnetic penguin op- 
erators^], for m s <C mi, 

Heff = ^V tb V t * s CM0 7 . (5) 
Here we use the notation of Grinstein et al. [12] 



7 = -^mbSa^F^^^b, 
16tt 2 m 2 

we denote it by a blob in Fig.l.a. 

The coefficient of 7 includes the QCD corrections, at the scale of W mass, 



(6) 



C 7 (m w ) 



x 



24(1 -xY 



[(1 - x)(8x 2 + 5x-7) + Qx(3x - 2)lnx], 



(7) 



where x = mf/m^, and when it runs down to low scale \x = rrib, it turns out to be 



C 7 {m b ) = y 



23 



C 7 (m w ) - j^(v 23 - 1) 



29 28 

(y 23 _ 1) 



189 



(8) 



where y = a s (mb)/a s (m w ). We have neglected the mixing of 7 with other operators 
which give small effects. 

Now we write down the amplitude of Fig.l.a as 



M n = 



dx\dy 



i<Pd;<Pb, 



C 



x -\Tr 



d -m Dt )fa fll/ (l + -f 5 )k v rj* 1 "^ - yj +m h ) la ^ + m Bc ) lbl a 



1 1 



+ Tr 



W " m Dt )ef iM - K„(l + TiW + m Bc h 5l a \ (9), 



1 1 



where rj is the polarization vector of the photon, x\ and y\ are the momentum fractions 
shared by charm in B c and D*, respectively. The denominators are 



D 3 



(1 - yi)(m 2 Bc - m^yi) - m 
(1 -xi)(m^. -m\x x ) 
(xi - yi)(xim| - m^yi), 



(10) 



and the constant C is 



2^2^ 



'4tt ^2 



11) 
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The result is 



M a = ^a^VV^/l™ 



(77* ■ e*p ■ k — rf • pe* ■ fyfy 



"fpeng 



(12) 



with the formf actors 

fT* = fiT 9 = C f 1 dx 1 dy 1 5(x 1 - ^-)5{ Vl - 1 + Cj3 .)4 

Jo mn,_ 



X 



m Bc (l - yi)(m Bc - 2m D .) - m b (2m Bc - m D *) 



1 1 



1 1 



D 1 D S ™B e m D .(l-x 1 ) — —\,(U) 



2.2 The weak annihilation contribution 



Being made up of two different heavy flavors, B c meson is also the unique place to probe 
the weak annihilation mechanism. 

Using the formalism developed by Cheng et al. [3], the amplitude of Fig.l.b are found 
to be 



M b = te, ua ^ V *»q a kPfr m + (V* • e*p • k — rf • pe* ■ k)f. 



anni 
2 ' 



(14) 



with 



fanni _ 
f 2 = -Kdi 



e s _e^\ m D * 



m D *m Bc 



ml 



m 



2 ' 
D* 



+ 



m D *m Bc 



2 ' 
D* 



(15) 

,m s m c j m Bc \m c m b ; 
where k = eGpV cb V*J B J B * / y/2, and the is the constituent quark mass. The parameter 
ai, corresponding to the one appearing in the nonleptonic B decays, includes the QCD 
corrections to the four Fermion operator (cb)v-A{sc)v-A- In nonleptonic B decays, a\ is 
extracted to be 1.01 from the recent CLEO data on B -> D^n(p) and B -> J/W[13] 
by Cheng et al.[3}. Here we take a\ = 1 in our numerical results. 



3 Numerical Results 

In order to have a numerical estimation we adopt the following parameters. 
1. Meson mass and the constituent quark mass 



Mo* = 2.11GeV, m b = A.7GeV, m c = 1.6GeV, m s = 510MeV 
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from the Particle Data Group[14]. mg c has been estimated to be about 6.27GeV in the 
literature [15, 16]. 

2. Decay constants of B c and D*. Up to now f Ds has been reported by three groups 

232 ± 45 ± 20 ± 48MeU, S.Aoki et a/., [17] 
f D: = { 344 ± 37 ± 52 ± 42MeU, D.Acosta et a/., [18] 
430 ±\§, ±40MeV, J.Z.Bai et a?., [19] 

where the first errors are statistical, the second are systematic, and the third are uncer- 
tainties involved in extracting the branching fraction B(D+ — > yU + z/ M ). From heavy quark 
symmetry, f D * = f Ds) we take f D * = 3AAMeV. 

The pseudoscalar decay constant fs c is related to the ground-state cb wave function at 
the origin by the Van Royen-Weisskopf formula[20] modified for color 



h 



f Bc = J— I *ioo(0) |= J— | flio(0) | • 
V m B c V nm B c 

Using the recent results for R 10 derived by Quigg[21]. One finds 

500MeU, (Buchmuller - Tye potential[22}) 
510MeU, (Power — law potential[23\) 
479 MeV, (Logarithmic potential[24\) 
696MeU, (Cornell potential[2h\) 

we will take fs c = 500MeU for our numerical estimation. 
3.CKM elements and strong couplings 

For quark mixing matrix elements, we will take V cb = 0.040 [26], | V ts \— V cb , \ V cs \— 0.9745 
and V t b = 1. The QCD coupling constant a s (/i) at any renormalization scale ji can be 
calculated from a s (M z ) = 0.117 via 

( v _ a s (M z ) 
^ l-^-^ln(M z /,) 

and = 11 - \n s . We get C 7 (M W ) = -0.1953 and C 7 (M b ) = -0.2939. 

With the values given above for various quantities, we can compute the form factors, 
the decay rates and the branching ratios. 

The results are 

flpeng = i m x ^-9^ jjpeng = ffpeng^ f anm = L?24 x ^-9^ f anm = ^373 x ^-W 



in units of GeV 1 . 

The ratios of long- and short distance contributions to the form factors /i and f 2 are 

fanni fanni 
Jl -, A nA J2 



npeng 



1.474, 



fipeng 
J2 



-0.710, 



Jl J2 

which are independent of fs c and f D *, and expected to be reliable. 
From the amplitude formula eq(12)(14), we get 



rm 



[m 



B,- 



m 



D' 



32nm 3 Bc 



ifi+f: 



T(B C -> D* sl ) = 



We finally obtain the decay rates 

4.68 x 10~ 18 GeV, only penguin 

6.25 x 10~ 18 GeV, only annihilation 
1.45 x 10~ 17 GeV, penguin + annihilation 

The lifetime of B c is predicted to be tb c = 1.35 ± .lbps by Quigg[27], but other theorists 
report smaller values tb c = OAps ~ 0.7ps[27]. We estimate the branching ratio Br(B c — > 
-D*7) as a function of r Bc - The results are given in Table 1. 



Br(B c ^D* sl ) 


t Bc = 0Aps 


t Bc = 0.7ps 


r Bc = l.Ops 


r Bc = l.SBps 


g r peng 


2.85xl0~ 6 


4.98xl0~ 6 


7.11xl0~ 6 


9.61xl0~ 6 


j^^anni 


3.80xl0~ 6 


6.65xl0- 6 


9.50xl0~ 6 


1.28xl0~ 5 


T3„.total 
DT 


8.83xl0~ 6 


1.54xl0~ 5 


2.21xl0~ 5 


2.98xl0~ 5 



Table. 1 



4 DISCUSSION AND CONCLUSION 

We have studied two mechanism contributing to the process B c — > D*^y. The short 
distance one(Fig.l.a) induced by electro-magnetic penguin has been estimated within 
PQCD framework which is suitable for the process involved large energy release and 
hard scattering. The distribution functions (pn( x ) coincide with the usual nonrelativistic 
Coulomb wave functions of cb and cs systems, and expected to work well in the present 
case. 

An unique nature of B c decays is that "the spectator" is heavy. Assume that the 
velocities of the heavy spectator are equal to Vi of B c (before scattering ) and to v 2 of 
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the final meson(after scattering ). In two body B c decays, v\ • v<i is always large. For 
B c — ■> -D*7, v± • v 2 — 1.7. Certainly, there is no phase-space for the propagators appearing 
in Fig.l.a to be on-shell, so the imaginary part of M a is absent. This is different from 
the situations in ref . [8] . Here we recall that the momentum square of the hard scattering 
exchanged by gluon is about 3.6GeV 2 , which is large enough for PQCD analysis. The hard 
scattering process cannot be included conveniently in the soft hadronic process described 
by the wavefunction of the final bound state, that is an important reason why we cannot 
apply the commonly used models with spectator ansatz, for example, HQET[7], BSW, 
ISGW[28] models, to the two body B c decays. That is also one of the reasons why the 
commonly used models cannot reliably predict the processes of B — > two light hadrons 
or 7+light hadron, since the large mass scale of the decaying B meson and the nearly 
massless final states imply that the transition amplitude is governed by hard process. 
Hence we conclude that the spectator ansatz is broken down in such kind of process and 
it is misleading to apply straightforwardly the commonly used models with spectator 
ansatz to such processes discussed above. Another competing mechanism is the weak 
annihilation. We find that it is equally important as the former ones. The situation is 
different from that of the radiative weak B ± decays which are overwhelingly dominated 
by electromagnetic penguin contributions. The results stem from at least two reasons: 

1. The compact size of B c meson enhances the importance of the annihilation decays. 

2. In B c — > D*7 case, the CKM factor of weak annihilation contribution is | V c f,V cs |, but 
in B ± — > K*^( case, the CKM part is | V u i,V us | which is much smaller than the former. 

We have neglected the contribution of the vector meson dominance(VMD). Using 
the methods developed in [29], combined with the PQCD calculation of the form fac- 
tors of the process B c — > D*J/^{^') (following the conversion J/^>(^') — > 7 ), we find 
that the VMD contribution is negligibly small, due to the small J/^/(^') — > 7 coupling 
e/#j/tf(*')7 ^0.025(0.016). 

Finally, we want to say a few words about the possibility of observing the interesting 
process at Tevatron and at the CERN Large Hadron Collider(LHC). The numbers of B c 
produced at Tevatron and LHC are estimated to be [30] 16000 (for 25Pb~ 1 integrated 
luminosities with cuts of P T (B C ) > lOGeV, \y(B c )\ < 1 ) and 2.1x10 s (for 100/& 1 inte- 
grated luminosities with cuts of Pt(B c ) > 20GeV, \y(B c )\ < 2.5 ), respectively. Using the 
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numbers and Br(B c — > D*j), we find that the decay channel is unobservable at Tevatron, 
but more than one thousand events of interest will be produced at LHC, so it can be well 
studied at LHC in the near future. 
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Table and Figure Captions 



Table. 1. Br(B c ^D*'~f) is displayed for varous B c lifetime. The second raw is for the 
contribution of electromagnetic penguin. The third raw is for the contribution of an- 
nihilation. The final line is for the total results. Fig.l.a. The leading penguin graphs 
contribute to B c — >-D*7- 

Fig.l.b. The leading W-annihilations to B c — >-D*7- Contributions due to photon emis- 
sion from other quarks are denoted by ellipses. 
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(b) 



Fig.l 
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